On the absence of the usual weak-field limit, and the impossibility of embedding some 
known solutions for isolated masses in cosmologies with f(R) dark energy 
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The problem of matching different regions of spacetime in order to construct inhomogeneous 
cosmological models is investigated in the context of Lagrangian theories of gravity constructed 
from general analytic functions f{R), and from non-analytic theories with f{R) — R". In all of 
the cases studied, we find that it is impossible to satisfy the required junction conditions without 
the large-scale behaviour reducing to that expected from Einstein's equations with a cosmological 
constant. For theories with analytic f{R) this suggests that the usual treatment of weak-field systems 
may not be compatible with late-time acceleration driven by anything other than a constant term 
of the form /(O), which acts like a cosmological constant. For theories with f{R) — R" we find 
that no known spherically symmetric vacuum solutions can be matched to an expanding FLRW 
background. This includes the absence of any Einstein-Straus-like embeddings of the Schwarzschild 
exterior solution in FLRW spacetimes. 
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I. INTRODUCTION 

Fourth-order theories of gravity have recently at- 
tracted a considerable amount of attention as they ad- 
mit Friedmann-Lemaitre-Robertson- Walker (FLRW) so- 
lutions that can accelerate at late times without the pres- 
ence of any exotic fluids. R then follows that the appar- 
ent need for dark energy could simply be due to an in- 
appropriate application of Einstein's equations to scales 
beyond those within which they have been thoroughly 
tested. This idea is compelling, but constitutes a rad- 
ical shift from the standard approach to cosmology. R 
must therefore be carefully studied in order to ensure 
that the consequences of deviating away from Einstein's 
theory are fully understood. In this paper we attempt to 
contribute to this understanding by considering the con- 
struction of inhomogeneous cosmological models within 
the framework of f{R) theories of gravity. 

The gravitational fields around isolated objects, and 
the FLRW solutions of f{R) theories, have both been ex- 
tensively studied in the literature (see [iHl] for reviews). 
Here we do not intend to contribute further to the study 
of cither of these fields, but instead to the ways in which 
one can construct cosmological models that contain mas- 
sive astrophysical bodies. This will be done by attempt- 
ing to match together existing solutions. In particular, 
we will attempt to construct 'Swiss cheese' models by 
matching spherically symmetric vacuum solutions with 
FLRW solutions, as well as constructing 'lattice models' 
by matching together large numbers of the spacetimes 
associated with regularly spaced astrophysical bodies. 

The motivation for this study is to understand both 
the effect of cosmological expansion on the gravitational 
fields of astrophysical bodies, as well as the large-scale 
expansion that emerges in a universe with large den- 
sity contrasts. These questions have been carefully stud- 



ied in Einstein's theory, where the above constructions 
have proven to be useful devices for understanding them. 
Fourth-order theories are considerably more complicated 
than Einstein's theory, but by applying the same con- 
structions we should expect to gain some insights into 
these questions. These extra complications include the 
absence of Birkhoff's theorem, so that spherically sym- 
metric vacuum spacetimes are not unique 0, as well as 
more complicated junction conditions [7|. 

Further motivation for this study comes specifically 
from the work of Mignemi and Wiltshire Q . These au- 
thors used a dynamical systems approach to perform a 
non-pcrturbative study of the static, spherically symmet- 
ric solutions of analytic f{R) theories. They found that 
these solutions are generically not asymptotically flat, 
and that boundary conditions could therefore be impor- 
tant in determining the gravitational fields of isolated 
massive bodies. Similar results have been found for non- 
analytic /(i?) theories @. These effects arc entirely ab- 
sent if one assumes asymptotic flatness from the begin- 
ning, as is standard in most approaches to studying weak 
gravitational fields [l^, El • The construction of inhomo- 
geneous cosmological models, as outlined above, provides 
a way to implement appropriate boundary conditions, 
and therefore allows the validity of standard weak-field 
approaches to be investigated. 

At present, much of the current literature assumes that 
in /(i?) theories the evolution of the FLRW 'background' 
cosmology proceeds independently of the growth of struc- 
ture within it. The motivation for this within Einstein's 
theory comes, in large part, from the studies of inhomo- 
geneous cosmologies, such as those discussed above. R 
also comes, however, from the correspondence between 
Newtonian cosmology and the FLRW solutions of Ein- 
stein's equations during dust domination: The rate at 
which nearby astrophysical bodies fall away from each 
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other can be considered as being due to a Newtonian 
force (up to the usual accuracy this imphes), or due to 
the expansion of the universe. Both are reasonable de- 
scriptions on small enough scales. If one attempts to use 
f{R) as an explanation of dark energy, however, then one 
wants the cosmological expansion to be different to that 
of a dust dominated universe. The usual interpretation 
of the motion of nearby astrophysical bodies as being de- 
scribable (up to some accuracy) within Newtonian theory 
is therefore lost, and the intuition we have gained on this 
subject from studying the solutions of Einstein's equa- 
tions must be re-evaluated. 



II. f{R) THEORIES OF GRAVITY 

The action for f{R) theories is given by replacing the 
Ricci scalar, R, in the Einstein-Hilbert action by a func- 
tion of the Ricci scalar, f{R), so that the gravitational 
Lagrangian density is £ = f{R). 

Including a matter action, and extremizing with re- 
spect to the metric, the field equations for these theories 
can then be written as 



r — 



T. 



(1) 



A thorough investigation of this problem is a very dif- 
ficult proposition, as in order to evaluate the existence or 
not of a weak-field limit, and the emergence of FLRW-like 
behaviour on large scales, one cannot begin by assuming 
the existence of either of these things. Any realistic inves- 
tigation, however, needs to make some assumptions, and 
here wc will begin by assuming that the gravitational 
fields around astrophysical bodies can be described by 
known solutions (either weak- field or exact). We will 
then proceed to see which FLRW solutions these can be 
matched with, or which FLRW behaviours emerge, given 
this assumption. This approach will not allow all of the 
questions posed above to be answered fully, but will al- 
low us to show that some situations that are possible in 
Einstein's theory are not possible in all f{R) theories. 



In Section|lT]we introduce f{R) theories of gravity, and 
in Section Uni we discuss what is meant by 'the weak-field 
limit'. Then in Section llVI we discuss the junction con- 
ditions that need to be satisfied when matching together 
different solutions in f{R) theories. In Section |V] we at- 
tempt to make a Swiss-cheese-like construction in which 
we match the usual weak-field solutions to FLRW solu- 
tions in theories with analytic f{R). Section IVTl contains 
an attempt within these same theories to match together 
many different weak-field regions to make a lattice model 
of the universe. In both Sections |V] and IVII it is found 
that these constructions are only possible if the large- 
scale behaviour is similar to that of solutions to Einstein's 
equations with a cosmological constant. In Section FVIII 
wc then proceed to try and match some known exact so- 
lutions, including here some theories with non-analytic 
f{R). We take three known exact solutions that describe 
spherically symmetric vacua in these theories, and try 
and match each to FLRW solutions. We find that in no 
situations can the junction conditions be satisfied at the 
boundaries between regions. Finally, in Section rVIIII we 
discuss our results. 



Throughout this paper we will us Greek letters to de- 
note spacetime indices, and Latin letters a, 6, c etc. to 
denote coordinates on a boundary. When it is required, 
the letters i, j, k etc. will be reserved for spatial indices. 



where 



/it/ 



/' 



(2) 

is the effective energy-momentum tensor of what we 
will call the "curvature fluid", and is the energy- 
momentum tensor of the standard matter fields. Primes 
denote differentiation with respect to R. 

That the field equations ([1]) are fourth-order in deriva- 
tives of the metric can be seen from the existence of the 
VyV^/' term in ([2]), a result which also follows directly 
from Lovelock's theorem. This is generally thought of as 
an undesirable feature in a Lagrangian based theory as 
it can lead to Ostrogradski instabilities in the solutions 
of the field equations. The f{R) theories, however, are a 
special case in which this instability can be avoided [l^ , 
due to the existence of an equivalence with scalar-tensor 
theories. In the special case f{R) = R\t can be seen that 
the fourth-order terms vanish, and Einstein's equations 
are recovered. 

The field equations ([T|) are automatically generally co- 
variant and Lorentz invariant as they are derived from a 
Lagrangian that is function of R only, and i? is a gener- 
ally covariant and locally Lorentz invariant scalar quan- 
tity. These same symmetries also guarantee that the left- 
hand side of ([T]) is covariantly conserved. The f{R) theo- 
ries therefore exhibit many of the key features of general 
relativity, while generalising Einstein's equations to allow 
new behaviour. This freedom has been shown to allow 
improved renormalisation of the gravitational interaction 
[l3| , as well as early universe inflation , and a possible 
explanation of dark energy. 



III. THE WEAK-FIELD LIMIT 

Before progressing further, let us begin by specifying 
exactly what we mean by the term 'weak- field limit'. We 
take this phrase to mean that in extended regions of the 
Universe that are small compared to the Hubble scale, 
but large compared to the Schwarzschild radius of any 
compact objects that may exist within it, that the geom- 
etry of spacetime within the region (but outside of the 
compact objects) can be well described by small fluctua- 
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tions around Minkowski space, such that 

— 'n^v ^ 1 (3) 

where rj^j_i, is the metric of Minkowski space, and there 
exists a coordinate system in which each of the compo- 
nents of h^ij is ^ 1 and slowly varying. The description 
given by Eq. ([3]), and the corresponding physics, is what 
is meant by 'the weak- field limit'. 

There are a number of points in this explanation that 
require further clarification. Firstly, what we mean by 
'Hubble scale' here is the quantity cH~^ when consid- 
ering space-like separations, and when considering 
time- like separations (here H is the Hubble constant, as 
measured by observers using the recessional velocity of 
nearby objects). For a region to be 'small' compared 
to the Hubble scale then means that any two points on 
the boundary of that region that are space-like separated 
should be ^ cH~^ apart, and that any two points that 
arc time-like separated should be ^ apart. This 

definition requires H to be reasonably uniform through- 
out each small region, which we will assume to be true. 
The criterion that these regions should be much larger 
than the Schwarzschild radius of any compact objects, 
and that Eq. ([3]) should not be taken to describe the 
regions inside (or near) compact objects, are simply in- 
tended to remove from our consideration the regions near 
black holes and neutron stars. 

Let us now further consider Eq. ([3]). The crucial point 
here is that the geometry of spacctime in the region un- 
der consideration can be taken to be close to that of 
Minkowski space. In this case one can decompose the ten- 
sor h^j^n according to how its various parts transform un- 
der spatial rotations in the background Minkowski space. 
In general, one can then write h^^j as (see [Tsj ) 

h^^dx^dx" = 2(l)c^df-2B,cdtdx''+2 {ip6ij + Hij) dx'dx\ 

The divergence of Bi and the trace of Hij can be set 
to zero by an appropriate choice of coordinates, and the 
divergence-less part of Bi and the trace-free part of Hij 
can be consistently ignored. This leaves 

ds^ ~ -(1 - 2(t))c^df + (1 + 2i})Sijdx'dx^ , (4) 

where and tp are both ^ 1 and slowly varying. We 
refer to this as 'the Newtonian limit' if (j) behaves like a 
Newtonian potential, and satisfies V^<^ ~ —AirGp. 

Finally, we can make the concepts of 'small' and 'slow' 
precise by introducing a dimensionless order-of-smallness 
parameter, e. Velocities, = dx"^ /dt, are then said to be 
'small' iiv/c^ 0{e), and quantities are said to be 'slowly 
varying' if acting on them with a time-derivative adds 
an extra 0{e) of smallness when compared to a spatial 
derivative (the order of smallness of time derivatives and 
velocities are expected to be the same as the evolution 
of gravitating systems are typically governed by the mo- 
tion of their constituents). From the field equations and 
equations of motion it can then be seen that the lowest 



order parts of (f> and iJj, and the energy density p, are 
given by 

The field equations and equations of motion within the 
region under consideration can then be expanded order- 
by-ordcr in e, with the 'weak field' limit of Eq. (|4]) cor- 
responding to the expansion up to O(e^). The ~ sign 
will be used in what follows to mean 'equal up to terms 
of 0(c') and smaller'. This is the same expansion in e 
that is routinely used in the standard parameterised post- 
Newtonian approach to gravitational physics in weak- 
field systems pjjf . 

One may note the differences between the perturba- 
tivc expansion outlined here, and the one that is rou- 
tinely used in cosmological perturbation theory about an 
FLRW background. This difference is intentional, and 
indeed necessary, for the study we are performing. Cos- 
mological perturbation theory is designed to be used on 
a variety of scales, all the way up to the scale of the 
cosmological horizon. This necessitates using a back- 
ground that is not Minkowski space and fields that are 
not slowly varying, as trying to model such a situation 
with a Minkowski background would involve 3-velocities 
that are of the same order of magnitude as the speed of 
light. Instead, here we are only interesting in modelling 
regions that are much smaller than the horizon size us- 
ing our perturbative framework, and then building up a 
cosmological model by joining together many of these re- 
gions. In this case the order-of-smallness parameter, e, 
can thought of (approximately) as the scale of the small 
region compared to the Hubble scale of the eventual cos- 
mological spacetime. In small regions such as this it is 
well known that relevant quantities become slowly vary- 
ing, and it is this slowness that is formally incorporated 
into the perturbative expansion in weak-field systems by 
allowing time-derivatives to add an extra power of e to a 
quantity. 

The weak-field limit of /(i?) theories of gravity 
has been studied extensively in the literature (see 
e.g. [H-Ill, and references therein), with the full post- 
Newtonian limit of theories with analytic /(i?) that ad- 
mit Minkowski space as a solution being found in [ll| . 
There the Lagrangian function is expanded in a Taylor 
series as 

/(i?) = /(O) + /'(0)i? + i/"(0)i?2 -f 0{R% (5) 

where primes denote differentiation with respect to i?. 
To the order required here the metric is then given by 
Eq. dH) with [m 

'/'=^(t/+^.ri?) , (6) 
'^-^(t^-^.f'i?) , (7) 
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where the functional dependence of /'(O) and /"(O) has 
been suppressed, and where U and R satisfy 
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and 



/' 



/ 



' 6/"- 



(8) 



(9) 



In these equations, and in what fohows, we have chosen 
to use geometrised units in which G = 1 = c. 

Assuming the existence of a weak-field limit, these the- 
ories can be seen to have a Newtonian limit if f"R <C U . 
Unlike in the PPN treatment, we will not insist that the 
solutions of Eqs. and ^ approach zero at asymptot- 
ically large distances, but will instead enforce boundary 
conditions using cosmology. 



IV. JUNCTION CONDITIONS 

In what follows we will be matching together different 
regions of spacetime in order to construct inhomogeneous 
cosmological models. This requires a set of junction con- 
ditions, analogous to the Israel junction conditions from 
general relativity [13], and is a problem that has been 
considered in f{R) theories of gravity by Deruelle, Sasaki 
and Sendouda Q. We will briefly recap the relevant re- 
sults from their work here, as it of central importance to 
our study. 

The central requirement in is that if one allows delta 
functions in the matter part of the field equations (i.e. if 
one allows matter fields to be localised on the boundary 
hyper-surface) , then delta functions should occur at most 
linearly in the parts of the field equations that involve 
geometry only. Here we are interested in the case in 
which there is no brane located at the boundary. We 
therefore require that there should be no delta functions 
in the part of the field equations containing geometry 
only. 

Now, in a Gaussian normal coordinate system, ds^ ~ 
dy'^ + jabdx'^dx'' , where the boundary is located at y = 0, 
the Ricci scalar can be written as 



R ^ 2dyK 



Kl^K*"" - -K' 



R, 



(10) 



where R is the Ricci curvature constructed from jab, the 
extrinsic curvature of the boundary is Kab = —^dy'^ab, 
and K and K*i^ are the trace and trace-free parts of this 
quantity, respectively. 

It can be seen from the field equations ([ij that R must 
be continuous at the boundary. This is because the cur- 
vature fluid contains terms like dyf'(R), which can be 
expanded as 



If R is not continuous then the second term above 
would contain a factor of S{y), this is not allowed un- 
less f"{R) ~ 0, which is just Einstein's equations. We 
can then see from Eq. pU)) that 'jab must also be contin- 
uous, otherwise Kab would contain a factor of S{y), and 
R would contain factors of {S{y))'^. This is not allowed, 
as Kab and R occur directly in the field equations. We 
therefore have that jab and R must both be continuous 
across the boundary. 

The yy and ya components of Eq. IT} are then given 

by 

dy [/' {Kab - Kjab) + labf'dyR] = 0. (12) 

Integrating this across the boundary one then finds 



[/' {Kab - Kjab) + labf"dyR]t = 0, 



(13) 



where the [. ..]!; notation means the difference of the 
quantity in the brackets on cither side of the bound- 
ary. Similarly, one can integrate R across the boundary 
to find, from Eq. pO|) that [R\t ~ 0, and hence that 

[2dyK - K2t,K*'''']t = 
of Eq. are then given by 



K^bK*"-"]! = 0. The trace and trace-free parts 



r{R) [dyR]t = 
f'{R) IKbt - 
[K]t = 

which, together with 

bab]t = , 

[R]t = , 



(14) 
(15) 
(16) 

(17) 
(18) 



dyf{R) = f"{R)dyR. 



(11) 



form the junction conditions in f{R) theories in which 
/''(i?) 7^ 0. For further details the reader is referred to 

a. 



V. MATCHING TO FLRW 

One of the oldest ways of trying to construct inhomo- 
geneous cosmological models is to join FLRW solutions, 
at some appropriate boundary, to the spherically sym- 
metric spacetimes that are expected to exist around in- 
dividual isolated objects. This was famously achieved by 
Einstein and Straus for the case of the Schwarzschild so- 
lution and the Einstein-dc Sitter universe [l^l- It is also 
possible to join the Lcmaitre-Tolman-Bondi solutions of 
Einstein's equations to FLRW at a spherical boundary 
[igt - These models are often referred to as 'Swiss cheese', 
as this is what the global structure starts to look like if 
one can keep removing regions of the FLRW 'cheese', 
and replacing it with either Schwarzschild or Lemaitre- 
Tolman-Bondi holes. The simplicity of this approach, 
and the degree to which it has influenced the develop- 
ment of inhomogeneous cosmology in general relativity. 
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makes it a natural place to begin studying the relation- 
ship between weak-field systems and cosmology in f{R) 
theories of gravity. 

Here we will consider an FLRW geometry given by 

dp 



ds^ = -dt^ + a\t) 



1 — A;f2 



'-de' 



sin^ 



(19) 

and that is filled with a perfect fiuid. Within this space- 
time we will excise a region interior to the sphere f = E, 
and replace it with a region of spacetime that is spheri- 
cally symmetric, and that is well described by the weak- 
field geometry given in Eq. Q. In this case it is conve- 
nient to write the spatial metric in spherical polar coor- 
dinates, so that Sijdx'-dx^ = dr^ -I- r'^d9'^ + sin^ (P-d(\? . 
We can then identify the angular coordinates in both re- 
gions, which we will refer to as Region I and Region II, 
respectively. 

Without loss of generality, we consider the boundary 
to be comoving with the fluid. As there are no spatial 
gradients in Region I, the boundary must be static with 
respect to the hypersurfaces of homogeneity that exist in 
the FLRW geometry. In Region II, however, the bound- 
ary is free to move in the radial direction. The first fun- 
damental form of the boundary, on either side, is then 
given by 



-/ildx^dx" 



-dt' 



- ( 1 - 20 - ) dt^ + (1 + 2i')j:'dn' 



where the boundary is at r = S in Region II, and where 
we have used the notation ~ to mean equal up to terms 
of post-Newtonian order (i.e. up to 0{e^)). The junction 
condition ([T7)) then gives the conditions 



(1 + tp)Y, ~ aE , 

dt ^2 



(20) 
(21) 



Now let us consider the extrinsic curvature. To calcu- 
late this we need to know the normal to the boundary, 
which is given in each region by 

(22) 



VI - fcf2 ' 
1, 

2" 



(23) 



The second fundamental form on the boundary is then 
given by 

dxf^ dx" 



(24) 



dx- dx" ' 

which for the two spacetimes we are considering is 
Ki,dx^dx' ^ , (^1 - g + 1(7 - i^R^ dn\ (25) 



Kiidx'^dx'' ~ ( jU^r + jpR,r - t ) dP 



(26) 



1 - —R + -T.' - —rR r dVt\ 
2f' Of I ' 



2/' 



where wc have already used the junction conditions ((20)) 
and (EH. 

The junction conditions (fTSl) and ([T6|) then give 



2U\, 



S ~ /'E ' 2/' E 



E2 ' /' 



f" Rr\ 



(27) 
(28) 



These look very much like the Friedmann equations de- 
rived from Einstein's equations, with the terms contain- 
ing the Newtonian potential U acting like the matter 
terms, and with the term involving the spatial curvature 
k playing its usual role. Here, however, we also have 
two additional terms containing derivatives of the Ricci 
scalar, R. These extra terms can be seen to contain all of 
the new behaviour that one obtains by generalising the 
gravitational Lagrangian from R to f{R)- 

So far we have only applied the junction conditions 
that exist in Einstein's equations: That the first and sec- 
ond fundamental forms on the boundary must be contin- 
uous if we are to avoid a surface layer of matter. Let us 
now apply the additional junction condition (|14p . The 
spacetime in Region I is homogeneous, so in this case we 
must have 



dyR = Rf- = 0. 



(29) 



Applying the junction condition (1141) then gives 

R,rk-0, (30) 

where we have used kT,'^ ^ 0(e2), as can be seen from 
Eq. ([27)) . This means that the last terms on the right- 
hand side of both Eq. (|?7)) and Eq. (|^5)) must vanish at 
0(e2), so that we are left with exactly the same equations 
as in Einstein's theory (up to the presence of /'(O) in 
the denominator of the terms involving U, which can be 
absorbed into constants, and the terms involving /(O) in 
Eqs. ^ and dS]), which act like A.). 

This treatment appears to show that the only Swiss 
cheese solutions that exist in f{R) theories of gravity 
must either have FLRW regions that behave in the same 
way they do in Einstein's theory (possibly with A, and 
up to possible small corrections), or it must be the case 
that the spacetime within the excised spheres cannot be 
described using the weak-field geometry given in Eq. (jj)) , 
and explained in Section IlIIl 

VI. MATCHING WITHOUT FLRW 

In the previous section we tried to join a region de- 
scribed by a weak-field geometry to a FLRW geometry. 
Here we take a different approach, and instead try and 
join together numerous different regions each of which is 
well described internally by a weak-field geometry of the 
form given in Eq. (j?)). Note that we do not require mul- 
tiple regions to be well described by the same weak-field 
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metric, but instead try and relate the different weak-field 
systems to each other by using the junction conditions 
discussed in Section ITVl 

We will proceed by taking a number of objects with 
the same mass and distributing them regularly in space. 
We will then take the domain of each object to consist 
of all the points in space that are closest to that object. 
What we mean by a 'regular distribution' here is then 
that the domains of all objects in the spacetime should be 
identical, up to translations and rotations. An example 
of this in i?^ is obtained by dividing the space up into 
a cubic lattice, and placing a mass at the center of each 
cube, the interior of which then acts as the domain of the 
mass at its center. 

These situations have been considered within the con- 
text of Einstein's equations by a number of authors. Ap- 
proximate matching schemes have been developed and 
studied for joining together Schwarzschild solutions in 
the spirit of the Wigner-Seitz construction from solid 
state physics 2^; 2 3l|. P erturbative treatments have also 
been attempted |24l42a | , as well as an exact treatment of 
the initial value proplem (27j , and numerical studies have 
been performed [28ll29| . This type of model seems ideally 
suited to a study of the relationship between weak-field 
systems and cosmology, as they allow one to perform a 
bottom-up construction of a cosmological model from the 
weak-field systems themselves. They also do not require 
the existence of FLRW geometry at the boundaries be- 



tween regions, and so one is allowed to move away from 
some of the restrictions of the Swiss cheese models dis- 
cussed in Section FVl 

From the symmetry of the situation the junction con- 
dition P7)) is automatically satisfied. We then need to 
consider the conditions ([TC)) and ([T51) . To do this wc need 
to know the extrinsic curvature of the boundary of each 
domain. If the unit normal to this hypersurface is given 
by — {nt;ni) then, in the geometry (j4]) this is given 
to lowest non-trivial order by j24{ 

K ^ijdx^ dx^ ~ {rit^t + nicp^i) dt^ -(- [rii^t + ^t,i) dx^dt 

+ (jiij — ^ipiUj + 5i,.j5^^ilj,uni) dx^dx-' . 

The magnitude of rit can be seen to be of 0{e) here, while 
rii is 0(1), as the condition u^n^ — 0, where is the 
4-vclocity of the boundary, gives rit = —UiU^, and is 

0(6). 

Wc must now transform A'^i^ into Kab using Eq. (|24p . 
To do this it is convenient to pick out a direction, which 
we call z, that is normal to the boundary at the point 
where it intersects a straight line that joins the positions 
of two neighbouring masses. The remaining two spatial 
directions, which we will use together with t as the in- 
trinsic coordinates on the boundary, will then be written 
using indices A,B,C etc.. The position of the bound- 
ary itself will be given by z = Z {t, x ^). The extrinsic 
curvature of the boundary is then [2J] 



Kabdx°'dx ~ —71,2 



(Z - cj)^, + Z^Acb^A) dt" + + {Z)\a) dx^dt 



+ [Z\AB - (VA. - <5^^Z|cV'|D)('5Ai3 + Z\AZ\B))dx^dx'' 



(31) 



where ' = u^^d^ = dt + Z^tdz and \a = m'^d^ = 
Oa + Z^aQz- Here is a space- like unit vector in the 
boundary, and we have taken u"^ = = to*. All quanti- 
ties in Eq. should be taken to be evaluated at the 
boundary. 

It can now be seen that in order to simultaneously sat- 
isfy Eqs. (fT5|) and ((T6|) . as well the reficctive symmetry 
about the boundary that is required by our regular dis- 
tribution of masses, we must have Kab = 0. From Eq. 
this can be seen to correspond to 24l |: 



Z 



(32) 



^1 



\AB 



{5AB + Z\AZ\B){n-yi'), (33) 



as well as (2'|^)' ~ {^)\a ~ ^' These equations govern 
the motion and shape of the boundary in these highly 
symmetric configurations. 



Finally, let us apply the junction conditions (jl4p and 
P^ . The latter of these is again automatically satisfied 
from the symmetry about the boundary. The former, 
however, gives n • Vi? = 0, which means that Eqs. ([5^ 
and (|33p become 



Z 



Z\AB ~ {^AB + Z\aZ\b)Z. (35) 

These are exactly the same expressions that are found us- 
ing Einstein's equations (up to the factors of /'~^, which 
can again be absorbed into constants, and terms involv- 
ing /(O), which again act like A). The solutions to Eqs. 
([M]) and ([55)1 are known to be either decelerating, or cor- 
respond to Minkowski space unless /(O) ^ 0. 

Once again, we find that for the weak field regions to 
exist we must have cosmological behaviour that is the 
same as in Einstein's theory (up to possible small correc- 
tions). Therefore, within this approach, the large-scale 
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behaviour must reduce to that expected from Einstein's 
equations, or one must rehnquish the usual weak-field 
description around astrophysical objects. 



VII. MATCHING EXACT SOLUTIONS 

We have so far considered joining weak-field geometries 
to either FLRW, or to each other, in theories in which 
f{R) is an analytic function. This has shown that accel- 
eration in the resulting cosmological model cannot occur 
in any new ways if the junction conditions given in Sec- 
tion HV] are to be obeyed. One must therefore either allow 
for gravitational fields to be rapidly varying, or give up 
on a description of the regions around astrophysical ob- 
jects as small fluctuations about Minkowski space. The 
latter of these two possibilities suggests that it may be 
useful for us to consider exact solutions. 

Unfortunately, the complexity of the field equations 
([T]) make exact solutions difficult to find. However we 
know that for all functions f{R) which are of class 
at R = and /(O) — while /'(O) 7^ 0, the Schwarzschild 
solution is the only vacuum solution with vanishing Ricci 
scalar [30| . It therefore seems natural to try and match a 
spherical region with Schwarzschild geometry to an exte- 
rior FLRW spacetime. In the context of Einstein's theory 
this corresponds to the well known Einstein-Straus ap- 
proach described earlier [l8| . Furthermore, if we restrict 
our considerations to f{R) = R^^^ then there are two 
known exact solutions (other than the vacuum solutions 
of Einstein's equations, that is, which are also solutions 
of these theories). A static spherically symmetric vac- 
uum solution with non-trivial asymptotics was found in 
, and a time-dependent spherically symmetric vacuum 
solution was found in In what follows, we will also 
try and join these two solutions to FLRW geometries. 



A. An Einstein-Straus-like Construction 

The constructions we consider here consist of point-like 
masses at the centre of otherwise empty spherical regions, 
whose geometry is described by the Schwarzschild metric, 
and that are embedded in FLRW geometry at appropri- 
ate boundaries. Such constructions were originally con- 
sidered by Einstein and Straus [l^ , and were introduced 
to address the question of whether or not the expansion 
of the universe can affect local mechanical phenomena, 
such as planetary orbits. Since the spacetime near the 
central mass is Schwarzschild, the planetary orbits are 
given by the usual time-like geodesies of this geometry, 
and the cosmic expansion does not affect them. Let us 
now investigate whether such a construction can be per- 
formed in f(R) gravity. 

We begin by writing the Schwarzschild solution as 

dr'^ 

ds"^ = ^A{r)dt^ + -j^ + r^ {d9^ + sin^ Odcj)^) , (36) 



where 



Air) = 1 - 



2M 



(37) 



Let us now try and embed this solution in an FLRW 
geometry, as specified in Eq. ([T9| . To do this, consider a 
boundary at f = S in the FLRW spacetime and r = S in 
the Schwarzschild solution. The first fundamental form 
on the boundary is then given in the vacuum region by 



-fabdx^dx^ 



A 



dt^ + T^d^f (38) 



and in the FLRW region by 

-fabdx^dx'' = -dp + a^ii)±^dn^, 



(39) 



where we have identified angular coordinates in the two 
different regions at the boundary and where dil^ = dd^ + 
sin^ 6 d(jp . The junction condition ([TT]) then gives 



S = a(t)S 



(40) 



(41) 



To calculate the second fundamental form we need the 
space-like unit vector normal to the boundary. In the 
vacuum region this is given by 




VA^ - S2 
while in the FLRW region it is 



S, 1,0,0 



(42) 



(43) 



The second fundamental form on the FLRW side of the 
boundary is then 



(44) 



Kabdx'^dx'' = a(£)SV 1 - kt^dn^, 
while on the vacuum side of the boundary it is given by 
3.^./! — .;4"^.i4 2./1 S 



2VA(A2 - £2)3/2 



' £2^3 
(A2 - S2) 



(45) 
(46) 



where all quantities should be evaluated at the boundary. 
Matching K^n at the boundary we obtain 



A' 



1 



A 



(l-A:E2/a2) 



(47) 
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Writing the above equation in the coordinates {i, f, 9, 4>), 
and using Eqns. (PU)) . (jH]) together with the form of the 
function A{r), we find 



2M 



(48) 



The left-hand side of the above equation is the usual def- 
inition of the Cahil-Macvitte function in FLRW space- 
times. 

Differentiating Eqn. pS]) with respect to t gives G^i — 
0. This implies that the total pressure (standard matter 
and curvature fluid) must vanish on the boundary, but 
as the pressure in the FLRW region is homogeneous, this 
means that the total pressure should vanish throughout 
the FLRW region. In this case, equating the time com- 
ponent of the extrinsic curvature will not give any new 
information. 

If we now impose the requirement that R should be 
the same on cither side of the boundary, from Eq. ([T8|. 
then we must have 



1 



1 d'^a{i) 



da{t) 
di 



. (49) 



The above equation combined with the condition of van- 
ishing total pressure, then implies vanishing total density 
(curvature fluid and standard matter) in the FLRW re- 
gion. Whatismore, putting i? = in equation ([!]) shows 
that the effective energy-momentum tensor of the curva- 
ture fluid must be proportional to g^i,. It then follows 
that the energy-momentum tensor of standard matter, 
T^, must also be proportional to g^j^, and so can only 
be a vacuum energy. It also follows that the FLRW region 
can only be Minkowski spacetime (in Milne coordinates, 
if A: = —1). Finally, from Eq. (fT4|) we see that the normal 
gradients automatically match identically, as i? = on 
both sides. 

It is curious that the Schwarzschild solution cannot be 
embedded in any FLRW spacetime (other than the triv- 
ial case of Minkowski space) in f{R) theories of gravity, 
unless the theory is linear in R. 



B. A Static Solution in 7?" Gravity 

An exact static, spherically symmetric vacuum solu- 
tion of f{R) = is given by (§] 

= -A{r)dt^ + + r2 (d6l2 + sin^ Odcf)'^) , (50) 
B(r) ^ ' 

where 

2.5(1 + 2^) Cl 

A{r) = r (1-^) + , 

f (1-5) 



B(r) 



Cl 



(1 - -1-4^2) (1 _ 2<5_ 2,52) I 



The Ricci scalar for this solution is 

65(1 + (5) 



i? = -- 



(1 - 25 - 252)a2j]2 



(51) 



We will now try and embed this solution in an FLRW 
geometry. To do this, consider a boundary at r = S in 
the vacuum region. The first fundamental form on the 
boundary is then given in the vacuum region by 



-fabdx°-dx^ 



A 



B 



dt^ + T.^dn^. 



(52) 

(53) 
(54) 



In the vacuum region the spacelike unit vector normal to 
the boundary is given by 



Matching the first fundamental forms then gives 
S = a(i)E 




-2,1,0,0 , 



VAB - £2 

The second fundamental form of the vacuum side is 

2BA rt^ + AB - AB^A r - 2ABt 
lytt = - 



2VA{AB - S2)3/2 



K„o = 



' S2B2A 
{AB - £2) 



(55) 

(56) 
(57) 



where all quantities should be evaluated at the boundary. 
The junction conditions and are then satisfied 

if 



E2 = AB 

ii = 



1 



B 



(l-fcS2/a2)_ 

{A^rB + B^rA) B{2A^rB + B^rA) 

~ 2^ 2(l-A;E2/a2) ' 

Consistency of these equations requires 

{AB^f - A^fB) 



(1 - fcr2) 



0. 



(58) 
(59) 

(60) 



Substitution from Eq. ([501 shows that this can be 
achieved only if 5 = or —1/2. 

If we now impose the requirement that R should be 
the same on either side of the boundary, from Eq. (|18p . 
then wc get 



1 d^a{i) 



1 f da{i)\ k 

1^ V di J ^ 1^ 

5{l + 5) 
(l-25-252)a(£)2s2- 



(61) 
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This strongly constrains the allowed form of a{t). Finally, 
from Eq. p^ . we find that we must have 



ABR, 



Vab - s2 



= 0, 



(62) 



as there are no spatial gradients in the FLRW region. 
This means that we must also require R^r = at the 
boundary in the vacuum region. This is only satisfied if 
(5 = or —1, as can be seen from the right-hand side of 

Eq. dnn). 

It is therefore the case that the junction conditions 
from Section IIVI can only be satisfied if (5 = 0, in which 
case the field equations ^ simply reduce to Einstein's 
equations. In this case the vacuum solution given in Eq. 
([50| reduces to the Schwarzschild solution, and Eq. ([6T|) 
no longer needs to be satisfied as f"~0, and the right- 
hand side of Eq. (fTTj) vanishes automatically. The vac- 
uum solution ([50)1 cannot, therefore, be used to model the 
gravitational field of an astrophysical object embedded in 
an FLRW universe in these theories, unless f{R) is lin- 
ear in R. This is despite the fact that this solution is the 
asymptotic attractor of all spherically symmetric, static, 
vacuum solutions of theories with /(i?) = R^~^^ Q, sug- 
gesting that the spacetime around astrophysical objects 
that are embedded in FLRW should be time-dependent. 



C. A Non-Static Solution in R" Gravity 

An exact solution for time-dependent, spherically sym- 
metric vacuum situations in f{R) = theories is 
given by Q: 

ds^ = -A{r)df + q^{t)B{r) [dr^ + r^dn^) , (63) 



where q{t) ~ t (i-*) , and 



A{r) 



1 _ £2 



C2 



2/a 



1(7+2(5-1 



B{r) = (^1 + -1 j A 

where <t^ = 1 — 26 + A5^ . The Ricci scalar in this case is 
given by 



R 



6(5(1 + (5)(1 + 25)(1- 45) 
(1 - d)H''A 



(64) 



Again, to match this solution with a FLRW exterior, con- 
sider a boundary at r = E in this solution. The first 
fundamental form on the boundary is then given in the 
vacuum region by 



jatdx^dx'' 



A - q^Bt^ ) dt^ + q^BY:^dVL^ 



Matching the first fundamental forms then gives 



l-^A-q^B^^ 



(65) 



(66) 
(67) 



and the unit vectors tangent and normal to the boundary 
are given by 



1 




(1,E,0,0), (68) 
(-E, 1,0,0). (69) 



A- Bq^t^ 



Calculating the second fundamental form for the match- 
ing surface we get 



2qq{B'^q^Y.^ -2ABY.) + q^Y?{2ArB-BrA)-2q'^ABY.-AAr 

J^ff — , : ; , ( 70) 

2\/ABg(yl-Bg2E2)3/2 

qi:{ABrT, + 2AB + 2qtB'^qj:) 
J^-ee — , = ■ (71) 



2VAB^JA- Bq^^^ 
Matching the Ricci scalar then gives 



1 d^a{i) 1 fda{i)\ k _ 5{1 + S){1 + 2S){1 ~ iS) 

a{i) dp ^(1)2 [Hrj (l-5)2tM(E) ■ 



Finally, the boundary condition n ■ Vi? = 0, gives 



4(i-a) 

2C2 / 1 — 1* \ " (1-3^-4^^) 



(72) 



S = ^{^ ^1 t' (I-*-) (73) 

aE2 (1-^)^(1 + ^)=^ Vl + f / 
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unless 6 = 1/4:, 0, —1 or —1/2, in which case n ■ Vi? = automatically. We can now construct an algebraic constraint 
for E by equating Ktt on either side of the boundary and using Eqs. (|66p and (|73p to remove a{t) and S. This gives 



(l]2+c2)(l-5)fc-2c2Cr2E 



- \/l 



(1 - 5)\Ja^ - 4cit'''u-^r 'l]8(l-5)/<T(5^2 _ ^2)2^2(1-5) 



(74) 



This equation must be satisfied at all times, but is clearly 
very difficult to solve for S directly. We can, however, 
perform a series expansion in C2. To zero order we then 
have the constraint 

V^l -fcS2 = l + 0(c2), (75) 

so that fc ~ 0. This says that the FLRW geometry that 
we are embedding within must be close to spatially flat. 
Using this in the first order equation then gives 

^-^C2=0 + O{cl), (76) 

so that the only possible solutions would appear to re- 
quire either tr = + 0(c2), or C2 = 0. The first of these 
possibilities requires S to be complex, which we are not 
interested in here, and the second is the requirement that 
the central mass vanishes. The matching of this latter sit- 
uation to FLRW is trivial, as the geometry in Eq. ([55)1 
can itself be seen to reduce to FLRW when C2 — s- 0. Once 
again we therefore appear to be unable to match solu- 
tions to FLRW, except when 5 = 0, or when the entire 
spacetime is FLRW anyway. 

The anomolous cases that remain are those in which 
6 = 1/4, —1 or —1/2, as in these cases Eq. ([73| can no 
longer be used. Of these 6 = —1 seems problematic as it 
corresponds to a Lagrangian density C =constant, which 
can hardly be said to be a Lagrangian for gravity at all. 
The cases 6 = —1/2 and 6=1/4 also seem problematic, 
as in these cases the field equations ([T]) contain terms 
that are ill-defined, with both numerator and denomina- 
tor reducing to zero. In all of these cases the Ricci scalar 
must vanish, so the only exterior FLRW geometry that 
one could match to would have to be Milne anyway. We 
do not, therefore, consider them to be of any interest for 
our current purposes. 

We therefore find that even for this non-trivial non- 
static solution, a matching with a FLRW exterior is not 
possible. This is true even though the solution itself ap- 
proaches FLRW asymptotically. 

VIII. DISCUSSION 

The idea that the late-time acceleration of the uni- 
verse could be explained by modifications of the Einstein- 
Hilbert action has recently attracted considerable inter- 
est, but a complete understanding of the consequences of 
such a radical shift away from the standard approach 



to cosmology is still far from complete. Here we at- 
tempt to add to this understanding by considering the 
construction of cosmological models that contain massive 
astrophysical bodies within the context of f{R) theories 
of gravity. Such constructions are key to understand- 
ing the effect cosmological expansion has on the gravita- 
tional fields of astrophysical bodies, as well as describing 
the large-scale expansion that emerges in a universe with 
large density contrasts. 

After a discussion of the junction conditions that need 
to be satisfied when matching together different solutions 
in /(i?) theories, a number of attempts were made to con- 
struct inhomogeneous cosmological models by matching 
different regions of spacetime. This was done both for 
theories with general analytic functions f{R) and for non- 
analytic theories with f{R) = -R". In all cases studied, 
it was found that it is impossible to satisfy the required 
junction conditions without the large-scale behaviour re- 
ducing to what is found from Einstein's equations with 
a cosmological constant. For theories with analytic f{R) 
this suggests that the usual treatment of weak-field sys- 
tems may not be compatible with late-time acceleration 
that is driven by anything other an effective cosmological 
constant given by /(O). For theories with f{R) = R" , we 
found that a number of well known spherically symmetric 
vacuum solutions could not be matched to an expanding 
FLRW background, including the well known Einstein- 
Straus-like embeddings of the Schwarzschild exterior so- 
lution in FLRW spacetimes. 

The absence of these constructions represents a cru- 
cial difference between f{R) theories and scalar-tensor 
theories of gravity. In the latter it is already known 
that Einstein-Straus-like embeddings are indeed possi- 
ble, both in cosmological and astro phy sical gravitational 
collapse scenarios (see for example [31|). This is true de- 
spite the extra junction conditions that are required in 
scalar-tensor theories, where the scalar field and its nor- 
mal derivative must be matched at the boundary. These 
two conditions may initially seem quite similar to the ex- 
tra conditions required in f{R) gravity {i.e. matching 
the Ricci scalar and its normal derivative). However, it 
turns out that the conditions in f{R) theories are much 
more restrictive, and give much stronger constraints on 
the spacetimes allowed on either side of the boundary. 
This is due to R taking a very specific form once an 
ansatz has been made for the metric (by specifying it 
should be give by Eq. or for example). 

These results are quite different to what is suggested 
by using linear perturbation theory around an FLRW 
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background in f{R) theories. In that case there seems to 
be httle impediment to including large density contrasts 
by allowing Sp to become large, while (j) and ip are re- 
quired to stay small. This difference could indicate that 
while the weak field solutions we have considered here are 
problematic, there may be some way of obtaining useful 
(approximate) spacetime geometries from the perturbed 
FLRW approach. However, it would seem somewhat odd 
to require knowledge about the geometry of the entire ob- 
servable universe in order to model the spacetime around 
a single astrophysical object (if this was the approach to 
be taken, it would also appear to mean that the cur- 
rent framework for interpreting precision tests of gravity 
would have to be re-thought). Alternatively, it may be 
the case that the difference between the bottom-up con- 
structions attempted here, and the top-down construc- 
tion of perturbed FLRW, could be indicating that cos- 
mological back- reaction is large in f{R) theories. This 
is certainly plausible, and should probably be expected 
when "screening mechanisms" such as the chameleon ef- 
fect come into play. 

As a final thought, it would be interesting to study 



the physical consequences of 'jumps' in the Ricci scalar 
and/or in the normal derivative of the Ricci scalar across 
the boundary. As is well known from the Israel junction 
conditions, a jump in the second fundamental form gives 
rise to surface stress-energy and surface tension on the 
matching surface that can, for example, be used to sta- 
bilize gravitational vacuum condensate stars ■ In a 
similar way, it is plausible that relaxing the extra match- 
ing conditions in f{R) theories could give rise to surface 
terms that might be of physical interest. This has been 
studied in the context of branc- world cosmology in [3| • 
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